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Abstract. We present new results of our X-ray mon- 
itoring of the transient binary supersoft X-ray source 
RX J0513. 9-6951 in the LMC and of our re-analysis of 
optical light curves obtained during the MACHO project. 
We have covered a complete X-ray outburst cycle with the 
ROSAT HRI detector. From the amplitude and timescale 
of the soft X-ray variability, tight limits are derived for the 
temporal behaviour of the white-dwarf radius and the ef- 
fective temperature of its envelope. A limit-cycle model is 
proposed to explain the observed optical and X-ray vari- 
ability, the characteristic timescales of the durations of 
the X-ray on and off states, and those of the transitions 
between both states. Our observations confirm that the 
radius changes of the white-dwarf envelope occur on the 
Kelvin-Helmholtz timescale. The duration of the X-ray on 
and off states is compatible with the viscous timescales of 
the inner and outer accretion disk, respectively. 
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1. Introduction 

Luminous supersoft X-ray sources (SSS) have been estab- 
lished as a new and distinct class of objects which are ob- 
servationally distinguished by their very soft X-ray spectra 
with temperatures on the order of 30 eV and luminosities 
of 10^^ — lO'^^ergs"^ (for recent reviews see Kahabka & 
van den Heuvel 1997; van Teesehng 1998). Several SSS 
have been identified as accreting close binaries with orbital 
periods of ~ 1 day or less. The most popular interpreta- 
tion of these systems involves a white dwarf which accretes 
matter via Roche-lobe overflow and an accretion disk at a 
rate of ~ 1 — 4 x 10~^MQ/yr, sufficient to permit stable 
quasi-steady nuclear shell-burning in the surface layers of 
the white dwarf, either because of thermal timescale mass 
transfer from a more massive (slightly evolved) main se- 
quence companion (van den Heuvel et al. 1992) or because 



of wind-driven mass transfer from a low-mass irradiated 
companion (van Teeseling & King 1998). 

2. The transient binary SSS RXJ0513.9 6951 

The luminous transient soft X-ray source RXJ0513.9- 
6951 (henceforth RX J0513) discovered in the ROSAT all- 
sky survey (Schaeidt et al. 1993) has been optically iden- 
tified as a high mass-transfer accreting binary system in 
the LMC (Cowley et al. 1993; PakuU et al. 1993) with an 
orbital period of 0.76 days (Crampton et al. 1996). Opti- 
cal monitoring has revealed that RX J0513 undergoes re- 
current low states at quasi-regular intervals, in which the 
optical brightness drops by 1 magnitude (Reinsch et al. 
1996; Southwell et al. 1996). The optical low-states are 
accompanied by a turn-on of the system in the soft X-ray 
range (Reinsch et al. 1996; Schaeidt 1996). 

The optical low states last for ^ 40 days and repeat 
about every 140-180 days. Such short time scales cannot 
be explained by the limit-cycle behaviour sketched by van 
den Heuvel et al. (1992) or by recurrent burning models 
(Fujimoto 1982). Within the framework of a shell-burning 
white dwarf an alternative explanation has been suggested 
by PakuU et al. (1993): The rather sudden changes in the 
soft X-ray fiux are the result of the direct response of the 
white dwarf to slight changes in the mass transfer rate. On 
the horizontal shell-burning branch, a small increase of the 
accretion rate may significantly affect the effective radius 
of the white dwarf envelope (Kato 1985). An increase of 
the photospheric radius by e.g. a factor of 4 implies that 
the effective temperature drops by a factor of 2. Given the 
extreme sensitivity of the ROSAT PSPC and HRI count 
rates on temperature, this in turn implies that the source 
may become undetectable although the bolometric lumi- 
nosity remains roughly the same (e.g. Heise et al. 1994). In 
this model of an expanding and contracting envelope the 
sudden drop of the optical flux, the colour variation, and 
the temporarily increased soft X-ray flux can be quantita- 
tively described by variations in the effective temperature 
of the hot central star and variations in the irradiation of 
the accretion disk (Reinsch et al. 1996). 
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Fig. 1. X-ray and optical variability of the transient 
binary RXJ0513. Upper panel: ROSAT HRI detected 
(filled circles) and upper limit count rates (triangles) of 
RX J0513 during a complete outburst cycle. Lower panel: 
light curves of different optical minima determined from 
CCD images taken for the MACHO project. Filled circles 
predicted optical light curves calculated from our X-ray 
data (see text). 

This model is supported by an independent estimate 
of the neutral hydrogen column density obtained with re- 
cent HST UV spectroscopy which constrains the X-ray 
luminosity during the on-state of RX J0513 to (2.5 — 9) x 
Iq37 ej-gg-i^ i somewhat below the Eddington limit, and 
confirms that the radius of the soft X-ray source is con- 
sistent with the radius of a non-expanded white dwarf 
(Gansicke et al. 1998). 

3. X-ray and optical monitoring 

In order to obtain tight limits for the temporal devel- 
opment of the radius and the effective temperature of 
the photosphere, we have monitored RXJ0513 with the 
ROSAT HRI detector at intervals of about two days cov- 
ering one complete X-ray outburst cycle (Fig. |l|, upper 
part). The source remained undetectable during the X- 
ray off state and showed a sudden increase of the soft 
X-ray flux by a factor of > 100 at the end of August 
1997, reaching maximum flux after ~ 5 days. The X-ray 
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Fig. 2. Absolute V magnitude as a function of the white 
dwarf photospheric radius Ri for the disk, the companion 
star, the white dwarf, and their sum. 



outburst lasted ~ 40 days and ended with a steep flux 
decline by again a factor of > 100 to non-detectability 
within 2 days. The slow flux variation during the X-ray 
outburst can be approximated by an exponential decline 
with a time constant of (34 ± 6) days. 

Detailed light curves of several optical low states 
are available from regular monitoring of RXJ0513 as a 
serendipitous source on CCD images taken for the MA- 
CHO project (SouthweU et al. 1996). In Fig. (lower 
part) the light curves of different observed optical min- 
ima are shifted in time such that the onset of the optical 
decline occurs at At = 0. MACHO data obtained quasi- 
simultaneously with our HRI monitoring indicate that the 
steep onset of the X-ray outburst occurred within 1-2 days 
after the beginning of the optical decline (W. Sutherland, 
private communication). 

The shape of the optical low states repeats fairly well 
with an uncertainty of some 5 days in the duration of 
the low state. Within this uncertainty, its length coincides 
with the duration of the X-ray outburst. Both, the X- 
ray and the optical light curves show a similar small flux 
gradient before the final fast transition to the X-ray off/ 
optical high state occurs. 

4. Light curve modelling 

To test whether the model of a contracting and expand- 
ing white dwarf can quantitatively explain both the X- 
ray light curve and the dips in the optical light curve we 
have calculated a predicted optical light curve from our 
X-ray data. First we used LTE white dwarf model atmo- 
sphere spectra (van Teeseling et al. 1994) to determine 
the photospheric radius as a function of the HRI count 
rate, where we assumed a distance of 50kpc, a bolomet- 
ric luminosity of lO^^ergs"^, and an absorption column 
of riH = 6 X lO^^cm-^ (Gansicke et al. 1998). Then we 
used the binary light curve code binary-] — h (van Teesel- 
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ing et al. 1998) to calculate the orbital average optical 
magnitude as a function of the photospheric radius Ri of 
the white dwarf. Since this code self-consistently calcu- 
lates the amount of irradiation from an extended white 
dwarf on the accretion disk and companion, including all 
possible shielding effects, this calculation is more accurate 
than the semi-analytic approach we used in Reinsch et al. 
(1996) and allows us to investigate how the results de- 
pend on the various parameters. We assume a mass ratio 
of M2/M1 = 2, an orbital separation a = 3.8 x 10"'^^ cm 
as appropriate for a quasi-main-sequence donor star and 
an orbital period P = 0.76 days, an orbital inclination of 
10°, a disk filling 80% of the average Roche-lobe radius, 
a uniform irradiation reprocessing efficiency of 77 = 0.5, a 
secondary temperature of 9000 K, and an accretion rate of 
3.4 X 10"^ Mq. 

Figure || shows the resulting total absolute V magni- 
tude as a function of and individual magnitudes of 
the disk, the companion star and the white dwarf. For 
Ri/a > 0.13 or i?i ^ 5 X 10^° cm, the expanded white 
dwarf is the dominant optical light source. With increas- 
ing Ri the disk first becomes brighter because of more 
effective irradiation, but becomes fainter again for i?i ^ 
2 X 10^° cm because an increasing part of the inner disk 
disappears inside the white dwarf envelope. 

In Figure 0, we have plotted the predicted optical light 
curve over the combined MACHO light curve. For data 
points with only an upper limit for the X-ray count rate, 
we assume a radius Ri ~ 4.5 x 10^*^ cm, which correctly 
reproduces the amplitude of the dip in the MACHO light 
curve. The 3(7 X-ray upper limit of 0.00014 cts/s for the 
X-ray off state requires a radius of i?i 1 x 10^° cm (us- 
ing L = 10^8 erg s-\ T <, 185000 K, Tih = 6 x lO^Ocm-^) 
during the optical bright state. Our calculations show that 
it is relatively easy to reproduce the observed optical dips 
from the X-ray data, with the correct amplitude and sur- 
prisingly accurate absolute magnitudes. It also illustrates 
that when the X-rays become detectable, the white dwarf 
photosphere has almost reached its minimum size and the 
optical light curve has almost reached the level of the faint 
phase plateau. The difference between the observed and 
predicted optical light curve immediately after the optical 
decline could be explained by the initial lack of an opti- 
cally thick inner disk after the white dwarf envelope has 
contracted to its minimal proportions. 



5. A limit-cycle model for RXJ0513 

The X-ray and optical light curves of RXJ0513 suggest 
that the system follows a kind of a limit cycle behaviour 
with four typical time-scales: the ~ 140 days of the optical 
high/ X-ray off state, the rapid transition (;^4 days) to 
the X-ray on/ optical low state, the ~ 30 days duration of 
this state, and, again, the rapid transition (^2 days) to 
the optical high/ X-ray off state. 



Limit-cycle model for RX J0513.9-6951 




(a) t = 0: white dwarf photosphere starts to contract, 

because accretion rate at core drops below critical value. 

(b) t = 5 days: WD photosphere reaches minimum radius, 

low-density inner disk starts to fill on viscous timescale 
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(c) t = 10 days: inner disk becomes optically thick for irradiation, 

accretion rate on WD rises. 

(d) 1 = 31 days: accretion rate on WD reaches critical value, 

photosphere starts to expand, accretion rate 
on WD core continues to rise. 

(e) t = 42 days: photosphere reaches maximum size, outer disk 

starts draining empty on viscous timescale because of 
very effective irradiation. 
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— (f) t ~ 160 days: WD photosphere starts to contract, 

because accretion rate at core drops below critical value. t 

Fig. 3. Limit-cycle model for RX J0513 (sec text). 



We propose here that this behaviour results from ex- 
pansion of the white dwarf photosphere in response to 
enhanced accretion onto the white dwarf, together with 
the reaction of the disk to increased irradiation by this 
expanded photosphere while the mass transfer from the 
companion star remains constant. At the (arbitrary) start 
of the cycle (Fig. ^(d)), let us assume we have an accre- 
tion disk supplying matter to a white dwarf with its non- 
expanded radius Ri ~ 10^ cm. Because the mass supply 
rate is close to the Eddington critical accretion rate Merit 
(Fujimoto 1982; Kato 1985), the white dwarf radius be- 
gins to expand (Fig. ^ (d)-(e)), as explained in Sect. 2 
above. This will, in turn, influence the disk temperature. 
An extended central source with radius i?i 3> H, where 
H is the scale height of the disk, produces a surface tem- 
perature Tirr at disk radius R in an optically thick disk 
(e.g. Adams et al. 1988) given by 
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Here p — Ri/ R, Ti is the temperature of the white dwarf 
photosphere, limb-darkening has been neglected, and 77 is 
the reprocessing efficiency of the disk surface. Increasing 



4 



K. Reinsch et al.: A limit-cycle model for the binary supersoft X-ray source RX J0513. 9-6951 



Ri has two effects: (i) at given radius R, the disk temper- 
ature rises approximately as T^^. (x Ri R~^ oc Ri , and (ii) 
the inner disk disappears in the hot envelope of the star 
(see also Sect. 4 above). 

The increase in disk temperature raises the mass- 
flow rate in the disk, since the disk viscosity coefficient 
V = acsH is increased. Here eg is the sound speed, and 
H = csi?'^/V(2GM)i/2 is the scale height of the disk. The 
disk is now no longer in a steady state, since the mass- flow 
rate within it exceeds the mass supply rate from the com- 
panion star at its outer edge. Its mass is therefore gradu- 
ally drained onto the white dwarf on a viscous timescale 
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where 10^^ cm is a characteristic disk radius (the radius 
of the Roche lobe is about 1.5 x lO" cm), H/R^ ~ 0.03, 
eg ~ 3 X lO^cms-i for Ti„ ~ 30 000K , and a ~ 0.1 is a 
typical value of the viscosity parameter. With the disk be- 
ing drained, the accretion rate onto the white dwarf even- 
tually drops below Merit, and the white dwarf reverts to its 
unexpanded state, the disk becomes cooler, and the sys- 
tem enters an optical low and X-ray on state (Fig. ^ (a)- 
(b)). The collapse of the expanded stellar envelope leaves 
the accretion disk with an inner hole of approximately the 
envelope radius which is gradually refilled by accretion 
from the outer disk. The disk temperature at a given ra- 
dius R decreases as T^.^ oc i?i but with i?i ~ 10^ cm 
in the optical low state the temperature at the edge of the 
hole (Eh =i 3 X 10^° cm) becomes ~ 32 000 K, similar to 
the outer disk temperature in the optical high state. The 
viscous time scale for refilling the hole, assuming again 
H/Rh ~ 0.03 and cs ~ 3 x 10^ cms'^ then is 
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This picture predicts a long X-ray off state and a shorter 
X-ray on state, with rapid (thermal-timescale) transitions 
between them, in quantitative agreement with what is ob- 
served. 



6. Short-term variability 

Besides the exponential decline, the soft X-ray flux of 
RXJ0513 shows significant variations of ~ ±0.1 HRI 
counts/s on timescales of hours to days (Fig. A time- 
series analysis of the detrended HRI count rates, however, 
reveals no clear periodicity in the range 0.1-10 days. The 
strongest signal is found at P = 0.4075 days but corre- 
sponds only to a 2a detection. We have arbitrarily phase- 
folded the residual fluxes on the suggested orbital period 
of ^ 0.76 days but find no obvious modulation of the light 
curve neither using the spectroscopic ephemeris (Cramp- 
ton et al. 1996) nor using the better defined photometric 
ephemeris (Alcock et al. 1996). 



Fig. 4. Observed variability of the soft X-ray flux of 
(3^XJ0513 during the X-ray on-state after subtraction of 
the exponential long-term flux decline shown in Fig. |l|. 

Although our analysis shows that variability on the 
orbital and possibly shorter timescales may be present, 
our data coverage is not sufficient to decide whether the 
flux variations are truly periodic or not. 

7. Conclusions 

The ROSAT HRI monitoring of a complete X-ray out- 
burst of the transient SSS RXJ0513 has shown that the 
transition between the X-ray on and off states occurs with 
a change of the soft X-ray flux by a factor of > 100 within 
2-4 days. In the model of an expanding and contracting 
white dwarf envelope this implies that the decrease of the 
effective radius by a factor of > 7 and the increase of the 
effective temperature by a factor of ^ 3 occur on the same 
time-scale during the X-ray turn-on and vice versa during 
the X-ray turn-off. 

The steepest intensity variations in the optical occur 

?i\a similar time-scale as the soft X-ray flux variations, 
his is consistent with our model that the optical variabil- 
ity is caused by the varying contribution of the accretion 
disk illumination by the expanding and contracting enve- 
lope of the white dwarf. 

The existence of typical time-scales of the optical high/ 
X-ray off state, of the optical low/ X-ray on state, and 
of the transition phases, and the fairly accurate repeti- 
tion of the optical lightcurve suggests that the observed 
variability is driven by limit-cycle behaviour. A possible 
self-maintained mechanism is the periodic change of the 
accretion disk viscosity in response to changes of the irra- 
diation by the hot central star. In this scenario, the mass- 
ffow rate at the surface of the white dwarf varies while the 
mass-transfer rate from the companion star remains con- 
stant. Our model can qualitatively explain the observed 
time-scales and requires no external mechanism like the 
episodic occurrence of star spots near the Li point to trig- 
ger the transition between the X-ray on and off states. 
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